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^ I Abstract: We study eigenfrequencies and propagator expansions for damped wave 

equations on compact manifolds. Under the assumption of geometric control, the 
propagator is shown to admit an expansion in terms of finitely many eigenmodes near 
the real axis, with an error term exponentially decaying in time. In the presence 
of a nondegenerate elliptic closed geodesic not meeting the support of the damping 
coefficient, we show that there exists a sequence of eigenfrequencies converging rapidly 
to the real axis. In the case of Zoll manifolds, we show that the propagator can be 
expanded in terms of the clusters of the eigenfrequencies in the entire spectral band. 
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^ : 1 Introduction and statement of results 

In this paper we shall study some problems arising in the spectral analysis of dissipative 
^T) '. wave equations on compact manifolds. In order to formulate the results, we shall begin 

by recalling some standard assumptions and hypotheses. 
J ; Let M be a compact connected smooth Riemannian manifold of dimension > 2, 

and let A be the corresponding Laplace-Beltrami operator. We consider the Cauchy 

problem for the wave equation with a damping term, 

.>'■ J (-A 2 " A + 2ia(x)D t )u = 0, (t J x)eR + xM J 

X 1 u\ t=0 = u E H\M), D t u\t=o = ut e L 2 (M). 1 ' } 

Here D t = dt/i, and the damping coefficient a is a bounded nonnegative function on 
M, which is not identically zero. For simplicity, we shall assume that a e C°°(M; R+). 
The motivation for studying this problem comes from the geometric control theory 
(see [1^] and j2]), where one is interested in the long time behaviour of solutions to 
(jl.lj) . in relation to the geometry of the underlying manifold and the control (damping) 
region. 

Associated with the evolution problem (jl.lj) is the solution operator U(t) = e**- 4 , 
t > 0, acting in the Hilbert space of the Cauchy data TC = H 1 x L 2 and mapping 
(uq, mi) G 7i to (u(t, •), D t u(t, •)). Here we equip TC with the norm 

u \ 1 1 V7 112 ill 1 |2 ,11 1 1 2 

= II VWo 1 1 £,2 + || U \\ l2 + || U\ 111,2. 
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The semigroup U{t) can be introduced by means of the Hille-Yosida theorem — see 
Section 2 for the corresponding resolvent estimates. The infinitesimal generator A is 
the operator 

A= ( —A 2ia(x) ) :H ^ H > 

with the domain D(A) = H 2 x H 1 . Here H s = H S (M) is the standard Sobolev space 
on M. It follows that the spectrum of A, Spec (A), is discrete. 
The energy of the solution u(x,t) of (jl.lj) at time t, 

E(u,t) = \ f (|Vu| 2 + \D t u\ 2 ) dx (1.2) 

is nonincreasing as t — > oo, and relations between the rate of the decay of the energy 
and the spectrum of A were studied in [TH] and [20]; below we shall recall some of 
these results. The asymptotic distribution of the eigenvalues of A has been studied 
by Sjostrand [22], and some further results in the case when the geodesic flow on M 
is periodic have been obtained in [T2]. We also remark that under the assumption of 
periodicity of the flow, much more precise results are possible in the two-dimensional 
analytic case, and we refer to ^1] , as well as to forthcoming papers in this series, for a 
detailed study of this problem. In the present work we shall be concerned with results, 
valid in all dimensions > 2. 

We notice that r £ C is an eigenvalue of A precisely when the equation 

(-A + 2ira(x) - r 2 ) u(x) = (1.3) 

has a nonvanishing smooth solution. The eigenvalues r will also be called the eigenfre- 
quencies. The multiplicity m(r ) e {1, 2 . . .} of an eigenfrequency r is defined as the 
rank of the spectral projection 

Kr = ^~ [ (T - A)-' dT } 
2717 J y 

where 7 is a sufficiently small circle centered at To- We shall let E TQ := IT^'H stand for 
the generalized eigenspace corresponding to To. 

It follows easily from (jl.3)) that if r is an eigenfrequency, then 

< Imr < 2|| a ||oo, 

and since (jl.3|) is invariant under the map (r, u) — > (— r, u), the set of the eigenfrequen- 
cies is symmetric with respect to the reflection in the imaginary axis. We remark that 
t = is an eigenfrequency corresponding to the constant solution of (|1.3[) . and in jTH] 
it was proved that r = is the only eigenfrequency with the vanishing imaginary part. 
Let us introduce a strongly continuous family of operators 

U(t) : L 2 -»■ H 1 , t> 0, (1.4) 
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which takes a vector / G L 2 to the first component of U (t) applied to a vector of the 
form (0,/). The mapping in (jl.4j) will be referred to as the propagator. A simple 
computation shows that in terms of U(t), the matrix representation of the semigroup 
U (i) is given by 

jjffy - ( v(t) U{t) \ 

U{t) ~ \ D t V(t) D t U(t) ) H ^ H > 

where V(t)f = D t U(t)f - 2iU(t)(af). 

In the self-adjoint case, i.e. when a = 0, it is true that U(t) = isint^—A/yf—A, 
and we have a Fourier expansion, 

sin t a/— A . ^-^ sinA ? t . . . , . l2 , , 
/— t- /O) = y. , <Pj(x), -Atpj{x)=X j tpj{x), 

V " A A ?6 Spcc(-A) ^ 

where the convergence is absolute in the case of smooth data. In this paper, we shall 
be interested in eigenfunction expansions of the propagator U(t), as t — > oo, in the 
case of the nonvanishing damping term. We are therefore concerned with relations 
between the long time behaviour of U(t) and the spectrum of A. This problem is 
closely related to that of completeness and summability of generalized eigenfunction 
expansions of the operator A. The classical theory of non- self adjoint operators, as 
described in [TT] and ^7j, shows that the space of finite linear combinations of the 
generalized eigenfunctions of A, © r6 s pec (^) E T , is dense in TC. However, as emphasized 
in [1], the lack of orthogonality in the direct sum above indicates that the generalized 
eigenfunctions may still fail to form a (Riesz) basis in TC. The question of basisness in 
the context of weak abstract non-selfadjoint perturbations of selfadjoint operators has 
been studied in |T7j , and when specialized to the case at hand, the results of [T7] apply 
in dimension one only, implying then the existence of a Riesz basis of the generalized 
eigenfunctions. See also [8 a and |21j . Here we shall be concerned with the case n > 2, 
and as well as using the classical non-selfadjoint theory, we shall also rely on direct 
microlocal methods. 

In addition to [T7], in this paper we have been motivated by recent studies of 
expansions for propagators on unbounded domains in terms of the resonance states — 
see j2H] for an overview of this work and further references. Our Theorems 1.2 and 1.3 
below complement the corresponding results in the theory of resonances in a situation, 
when the decay of states is caused by a direct dissipative mechanism on a compact 
manifold. In Theorem 1.4 we state a result which, to the best of our knowledge, has 
no analogue for problems on unbounded domains. 

When p(x, £) = £ 2 is the principal symbol of —A defined on T*M, we introduce 
the Hamilton vector field H p , and recall that the integral curves of the Hamilton flow 
exp (tHp) : p _1 (l) — > p _1 (l) are mapped via the projection tt : T*M — > M to geodesies 
on M. When T > 0, we put 

1 f T 

(cl)t = — / a(exp (tH p ) dt on p 1 (1). (1.5) 
T Jo 
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We define also 



0<A(T)= inf ;<o) Tl 



which is a bounded continuous function on R + . We recall from jTHI that the limit 

A(oo) : = lim A(T) (1.6) 

T — >00 

exists and it is true that A(oo) = sup T>0 A(T). 

In what follows we shall say that the geometric control condition holds if there exists 
a time T > such that any geodesic of length > T meets the open set {x; a(x) > 0}. 
In this case it has been established by J. Rauch and M. Taylor [20] that there is a 
uniform exponential decay of the energy E(u, t). A proof of this statement can also be 
extracted from the arguments given in Section 3 below. 

Our starting point is the following result, which is essentially due to [TE], and was 
also proved in 



Theorem 1.1 Assume that the geometric control condition holds. Then A(oo) > 0, 
and, for every e G (0,v4(oo)) ; there are at most finitely many eig en frequencies in the 
strip 

R + i[0,A(oo)-e]. (1.7) 

Our first result gives an expansion of the propagator U(t), as t — > oo, in terms of 
the eigenfrequencies in the strip (jl.7j) . provided that the geometric control condition is 
satisfied. This results is analogous to the corresponding one in scattering theory, giving 
an expansion of the wave group in the exterior of a nontrapping obstacle, in terms of 
exponentially decaying resonant modes. See [2E] arid references given there. 

Theorem 1.2 Assume that the geometric control condition holds. Let e G (0,A(oo)) 
be such that there are no eigenfrequencies with the imaginary part = A(oo) — e. Then 
there exists C = C(e) > such that 

U{t)= e itT p T (t) + E(t), t>0, 

reSpec(^l), Im-r<A(oo)-£ 

where p T (t) is a polynomial in t with values in C(L 2 ,H S ) for any s G R, and the 
remainder term satisfies 



We have 



E(t)\\c { L^)<Ce- (A( ^-^. {U 



ijf(x)dx 2 



p - oit)f = jMx)Tx> feL 



4 



In order to motivate the next result, let us recall that in [16 , G. Lebeau determined 
the optimal rate of the exponential decay of the energy E(u,t). See also [20j. Set 

D = infjlmr; r G Spec (A), |r| > 0}, 

and 

a = sup{/5 > 0; 3B > \/u G H, Vt > 0,E(u,t) < Be" pt E(u,0)}. 

Here the space TC of the Cauchy data is identified with solutions of (jl.lj) . It was proved 
in BS| that 

a = 2min(D,A(oo)). (1.9) 

In particular, it follows from (jl.9j) that if there exists a geodesic not intersecting the 
support of a, then a = A(oo) = 0. In PE], an explicit example of a surface of revolution 
was given, where there exists a closed geodesic not meeting supp(a), while D > 0, so 
that the eigenfrequencies are bounded away from the real axis. In the example the 
geodesic in question is hyperbolic. Our next result shows that if there exists a closed 
geodesic away from the damping region, which is elliptic, then we have a sequence of 
eigenfrequencies approaching rapidly the real axis. 

The precise assumptions on the geodesic will be formulated in terms of the associ- 
ated linearized Poincare map, and we shall pause here to recall this notion. Given a 
simple closed .Hp-trajectory 7 : [0,T] — > p _1 (l) (satisfying 7(0) = j(T)), the linearized 
Poincare map P 7 is defined as follows: let H C p _1 (l) be a smooth hypersurface in- 
tersecting 7 transversally at 7(0). Then P 7 is the differential at 7(0) of the smooth 
locally defined Poincare map: H — > H, obtained by following the flow of H p once 
around 7. Up to a symplectic conjugation, the mapping P 7 can be viewed as a linear 
symplectic transformation on R 2(n ~ 1 \ and we say that 7 is elliptic and nondegenerate 
if the spectrum of P 7 lies on the unit circle and 1 is not an eigenvalue. 

Given a nondegenerate elliptic closed .Hp-trajectory 7, with the eigenvalues of the 
corresponding P 7 of the form e ±lctj , j = 1, ... n — 1, and < (Xj < tc, we say that P 7 is 
iV-elementary if 

n-l 

kjCtj 7^ mod 2%Z 

3=1 

for each k = (k u . . . fc n _i) G Z™" 1 with < £™=i |%| < N. 

Theorem 1.3 Suppose that there exists a closed primitive H p -orbit 7 : [0,T] — ► p~ l (l), 
which is nondegenerate elliptic, and such that the associated mapping P 7 is N-elemen- 
tary for some N > 4. Assume that 7r(7([0,T])) D supp(a) = 0. Then there exists a 
sequence ijj), j = 1, 2 . . . , of eigenfrequencies such that |ReTj| — * 00, and < Im Tj < 
Ojv(l) |Re7j| _ , j — > oo ; for any N. 

Remark. This theorem can be viewed as an analogue of the results of [23) [23] . 
and j2S] m the theory of resonances. 
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Our last theorem is concerned with an expansion of the propagator U(t) on a Zoll 
manifold, i.e., a manifold all of whose geodesies are closed. We refer to [3] for the 
background and motivation for the study of such manifolds. In this case, the set of the 
eigenfrequencies displays a cluster structure, and the expansion holds in terms of the 
disjoint clusters in the entire spectral band. 

Theorem 1.4 Let M be a Zoll manifold, i.e. assume that exp (7rH p ) (x, £) = (x,£), 
(x, £) G Then 

1. There exists C > such that all eigenfrequencies r, with Rer > ; except for 
finitely many, are contained in the union of the rectangles 



a C , a C 

k + T - T' k + T + T 
4 k 4 k 



+ t[0,O(l)}, fc = l,2,... 



1.10) 



where a G Z is t/ie Maslov index of the closed H p -trajectory 



{exp(tH p )(x,0; t G [0,tt]}. 



= E I E e< V,*(t) + E e*X fc (f) J , t > 0. 

fe=l \rSSpec(.4),Te/fc TGSpcc(^), -re4 / 

, (L11) 

i/ere p Tj fc zs a polynomial in t with values in C(L 2 , H s ) for any s G R. TTie oirfer 
sitm m (jl.lljl converges absolutely in C(H 9 , L 2 ) for every 9 > and for each fixed 
t > 0. 



Remark. The first part of Theorem 1.4 has been established in |12j . This is a well- 
known and in general optimal result in the self-adjoint case, a = — see Chapter 29 

of ng. 

Remark. The rectangles in (|1.1UJI become disjoint only for /c large enough, k > k ^> 
1, and the expansion in (jl.llj) should be understood in the following sense: If we set 

Ukit)= E j tT Pr,k(t)+ E ^XfcW; fc>A: , 

TSSpec(^l), re/fe TESpec(A),-TeI k 

then the series 

oo 

E 

k=ko 

converges absolutely in C(H e , L 2 ), for any 9 > and each t > 0, and we have 
17(f) = E e^ r (t) + E^W iryC(H e ,L 2 ). 

reSpec(^l), |Rer|<fco+a/4 fc>fc 
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The plan of the paper is as follows. After some preliminaries in Section 2, Theorem 1.2 
is proved in Section 3. The proof relies on a simple conjugation by pseudodifferential 
operators, microlocally in the region of the phase space, where we lack ellipticity. We 
also adapt the classical argument of Morawetz JHJ to our situation, in order to obtain 
Theorem 1.2 in the case of the L 2 initial data. In Section 4 we study the case when 
the geodesic flow is periodic. Here we recall the methods of and show how they 
can be used to prove the second part of Theorem 1.4. Section 5 is devoted to the proof 
of Theorem 1.3. The proof depends on an a priori semiclassical resolvent estimate, 
and here we follow the idea of [17] and |22| of creating gaps in the spectrum of the 
unperturbed operator by means of an auxiliary finite rank perturbation. At this stage 
we also use the results and methods of P, and [23J. 

The results of this paper have been announced in [TB*] . 

Acknowledgements. I am very grateful to Maciej Zworski for numerous very helpful 
discussions on the subject of the present paper. I would also like to thank Nicolas Burq 
for pointing out the argument of Morawetz, a modification of which is used in Section 
3. The support of the Swedish Foundation for International Cooperation in Research 
and Higher Education (STINT), as well as of the MSRI postdoctoral fellowship, is 
gratefully acknowledged. 

2 Preliminaries 

In the beginning of this section we shall recall, for future reference, some resolvent 
estimates for A. In doing so we shall write D = {r G C; Imr G [0, 2| | a | |oo]}, and we 
recall that Spec (^4) C D. 

Proposition 2.1 There exists a constant C > such that 

Proof: Let us consider the following closed codimension one subspace of 7i = H 1 x L 2 , 

Ho = {(u , ux) G 7i] J 2ia(x)uo{x) dx — J ui(x) dx = 0}. 

Since a > is not identically zero, it is true that the inner product on 7i , 

((w ,Wi), {vq,Vi)) e := J Vu Vv dx + J u{v{dx 

defines a norm || • \\e there which is equivalent to the natural norm on 7i, || • ||^. 
When U = (u ,ui) G D(A), it is also true that AU = (ui, —Au + 2iaui) is in 7i , and 
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a simple computation shows that Hq contains © ^TeSpec(.4) E T as a (dense) subspace. 
When U E D(A) H Hq, we have 



< Im (AU,U) E <2\\a\\ 00 \\U 



|2 



so that when Imr < 0, we get Im ((A — r) U,U) E > |Imr| (U,U) E . It follows that 
there exists C > such that for any U ETCo, 

\\{r-A)- l U\\ n <-^—\\U\\ H , lmr<0. (2.2) 



llm 



T 



Now it is easily seen that the space E is of dimension one and is generated by (1,0). 
We have H = 7Y © E , also in the topological sense, and since 

||(r-^)-V|| w < — L-||y|| w , lmr<0, V E E , (2.3) 
|lmr| 

we conclude that (|2.1j) follows from (|2.2|) and (J2.3)) . in the case when Imr < 0. The 
proof is complete, for in the case when Imr > 2|| a H^, we argue in a similar way. □ 

Remark. A straightforward computation shows that 

(t-A)-i=( ^){2ia-T) -R{r)\ 
[T ^ \ R{r){2iaT - t 2 ) - 1 -R{r)r J ' 

where 

i?(r) = P(t)-\ P(r) = -A + 2iaT-T 2 . (2.5) 
Combining (|2.4j) with Proposition 2.1 we get 

1 1 R(r) \ \c(l 2 .l 2 ) < I I it , , Im r < 0. (2.6) 
|r] |lmr| 

and 

C 

II i2(r) IU(L2,L2) < — — — Imr > 2110 1100. (2.7) 

|r| (Imr -2\\a\\ 00 ) 

It follows from the proof of Proposition 2.1 that 

C 

I ( r — AY n | \r(H n) < t; — t — > r ^ D, n — 1,2, ... , 

and an application of the Hille-Yosida theorem shows that A generates a uniformly 
bounded strongly continuous semigroup U(t) = e ltA , t > 0, in TC, so that in particular 
U(t) E C(R7, ^(L 2 , iJ 1 )) n C 1 (RZjT, £(L 2 , L 2 )), in the strong sense. See 0. We have 
the following relation between U(t) and the resolvent of A, valid when x E D(A 2 ), 

roc— ia 

e i '- 4 x = x + — / e itr (r-^) _1 ^x— , t > 0, a > 0. (2.8) 
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The integral in ()2.8|) converges absolutely in TC. If we choose x = (0, /), / G H 2 , then 
(|2.4j) shows that the first component of (r — A) 1 Ax is equal to —rR{r)f. Thus, with 
the absolute convergence in L 2 we have, 

■ roo—iot 

U(t)f = — / e UT R(r)f dr, a > 0, / G tf 2 (M). (2.9) 

Notice that when / G L 2 (M), ()2.9j) still holds, with the integral in ()2.9j) converging in 
L 2 (Imr = -a;L 2 (M)). 

When establishing our results, it will be convenient to work in the semiclassical 
setting. We shall therefore perform a semiclassical reduction, following [22j and [T2*] . 
When considering the eigenvalue problem 

(-A + 2ia(x)r - t 2 )u(x) = 0, 

with Imr = 0(1) and Rer 3> 1, we write r = \fz/h, where < ft < 1, and z belongs 
to the fixed domain [a,/3] + i[— 7,7], 0<a<l</3<oo, 7>0. We are then led to 
study the ft-differential operator 

V = P + ihQ(z), (2.10) 

where P = —ft 2 A and Q(z) = 2a(x)yfz. We notice that P is essentially selfadjoint 
on C°°(M), with the domain of the selfadjoint realization being H 2 (M), and Q(z) is 
bounded and selfadjoint for z real positive. 

In the main part of the paper, we shall make use of some elementary semiclassical 
pseudodifferential calculus, and we digress here to recall some relevant notation. Let 
S m (T*M) = S^ (T*M) be the space of functions a(x,£; ft) on T*M x (0, ft ], ft > 0, 
which are smooth in (x, £ ) and such that 

^a(x,e;ft) = 0^(l)(l + |e|) m - |/31 , (x,0eT*M, ftG(0,ft ]. (2.11) 

When a depends on some additional parameters, we require the estimates (|2.11|) to 
hold uniformly with respect to these parameters. The formula for the classical ft- 
quantization, 

1 f f ( x —y,0 

a{x,hD x ]h)u{x) = j^-— J J e l ^^a(x, f; h)u(y) dyd£, 

then defines a class of ft-pseudodiffential operators on M, which will be denoted by 
Ph (S m )- 

For ft-dependent symbols, we write a G S™ if there exists a G S m , independent of 
ft, such that a — ao G hS™ 1 " 1 , and we say that a$ is the principal symbol of a and of 
the corresponding operator. When A G Op h (S' m ) has a principal symbol ao G S m , we 
write A = Op h (ao) and say that A is classical. 
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We shall finally recall the notion of the (semiclassical) wave front set, WFh(A), of 
an operator A G Op h (S m ). When p G T*M, we say that p g WF h (A), if the full 
symbol of A, for some choice of local coordinates near the projection of p, is of class 
S~°°' m := r\kh k S m near p. It follows from the definition that WFh(A) is a closed subset 
of T*M, WFh(AB) C WF h (A) n WF h (B), and WF h (A) = hr k A G Op h (S m ), for 
any k. In this paper, we shall work in a purely semiclassical setting, and consequently, 
only compact subsets of T*M will be important. The definition of WF^ given above 
will therefore be sufficient for our purposes. 



3 Proof of Theorem 1.2 

The main step in the proof of Theorem 1.2 consists of establishing the following result. 

Proposition 3.1 Assume that the Rauch-Taylor condition holds. Then for every e G 
(0, ^4(oo)) there exists C = C{e) > such that iflmr G [—1, A(oo) — e] and |Re t\ > C 
then 

IWH^^r^. (3.1) 

1 + T 



Proof: Performing the semiclassical reduction of Section 2, we can reformulate the 
statement of the proposition in the following form: for every e G (0, A(oo)) there exists 
C > and ho > such that for all h G (0, ho] and all 

z G tt(e) := {Re z G [a, /5],-2v / ReI< Imz/fc < 2v / Re^ (A(oo) -e)}, (3.2) 
it is true that 

IKV-z^Wc^^Kj. (3.3) 

Here P has been defined in (12. 10|l . The proof of ()3.3|) . sketched in ^Hj was a contradic- 
tion argument, and it relied upon a propagation estimate for a suitable semiclassical 
measure. Here we shall give a different proof, which has the merit of being more direct. 
In doing so, we shall use some arguments of [22] and [T2] . 

In what follows we shall assume, as we may, that Imz = O(h). Following the 
argument of Section 2 of [22], we conjugate V = P + ihQ(z) by an elliptic selfadjoint 
/i-pseudodifferential operator A = Op h (e 9 ), where g G S' _1 (T*M) is to be chosen. We 
have 

A' 1 (P + ihQ{z)) A = P + ihOp h (q(Rez) - H p g) + h 2 R(z). (3.4) 

Here q(z) = 2a^fz and R(z) G Op h (S°). Let us choose some 5 > 0, sufficiently 
small but fixed. Continuing to follow the argument of [22], we see that we can choose 
g = g(Rez) G S^ 1 such that 

q(Rez) — H p g = (g(Re z))t on ([a — 5, (5 + 5)) . 
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Here the trajectory average {q(Rez)) T is given by 

1 f T 

(q{Re z)) T = — J q(Re z) o exp (tH p ) dt, 

and T > is to be chosen later, large enough but fixed. We then consider the right 
hand side of (|3.4|) . 

V = V T = P + ihOp h (q T ) + h 2 R T (z), (3.5) 
where qr € S° is such that 

q T = {q(Rez)) T on - 5, (3 + 5]). 

We now claim that 

lim inf (q(Rez))T = 2y/RezA(oo), 

T^oop~ 1 (E) 

locally uniformly in E > 0. Indeed, we observe that the homogeneity property of the 
ifp-flow implies that 

inf (a) T = inf (a)^ T , 

P \E) P 

and choosing T large enough but fixed, depending on e > 0, we can therefore arrange 
that for q T in (|3.5j) . we have 



q T 



> 2VRei (A{oo) - e/2) on p _1 ([a - 5, /? + 5]). 



Take now a cutoff function < x ^ Co°(P _1 (( Q! ~~ ^/^ + 5))), such that x — 1 near 
p _1 ([a:, /?]). An application of the Garding inequality allows us to conclude that for 

z e n(e), 

(im (V - z) X u\ X u) > r]h\\ X u\\ 2 ~ 0(h°°) \\u\\ 2 , (3.6) 

where rj = rj(e) > is fixed and h > is small enough, depending on e. Here we let \ 
also stand for the corresponding quantization, and the inner product and the norm in 
(J3.6j) are taken in L 2 (M). On the other hand, we have 

(im (V - z)xu\xuj = Im ({x{V - z)u\xu) + {[V, xM***)) , 
and the absolute value of this expression does not exceed 

0{1)\\(V - z)u\\\\ X u\\ + 0(l)\\[V,x}u\\\\xu\\. (3.7) 

Let now xi £ C^(p~ 1 ((a — 5,(3 + 5)) be such that Xi — 1 near a ,P\) an d the 

support of Xi is contained in the interior of the set where x — 1- Then 

llfcxHI < ||[P,x](i-xiHI + 0(/niMI< owlKi-xiVH + ocoihll- 



n 



Now the operator V — z is elliptic away from p 1 ([a, /?]), and the semiclassical elliptic 
regularity shows that 

II (1 - Xi)u\\ < 0(1)|| (P - + 0(011 « II- (3-8) 

We get 

1 1 [V, X ]u\\ < 0{h) \\{V-z)u\\ + 0(h°°) \\u\\. (3.9) 
Combining (jH.fi j) together with (j3.7j) and (j3.9j) . we infer that 

||X«|| 2 < {V-z)u\\ WxuW + Oih^WuW 2 

< ^\\(V-z)u\\* + l\\ X u\nO(h°°)\\u\\ 2 , 

and therefore 

< ^j^ll (V-z)u\\ + 0(h°°)\\u\\. 
Combining this with (j3.8j) . with x in place of Xi, we get 

||w|| < ^t^II (V-z)u\\, ueH 2 {M), zeQ(e). 
This completes the proof if we take into account the uniform bounds 

A, A' 1 = 0,(1) :H S ^H S , s G R, 

where the if s -spaces have been equipped with their natural /i-dependent norms. □ 

It is now easy to obtain an expansion of U(t) in C(H e ,L 2 ), 9 > 0, as t — > oo. 
Indeed, if / G H 2 then 

R(r)(-A + 2 i ar-r 2 )f = f, 

which together with Proposition 3.1 implies that for | Re r | > Candlmr G [—l,A(oo) — 
e], we have 

\\R{r)\\c(H^)< {1 °^ )2 - (3.10) 
Interpolating between (j3.1j) and (j3.10j) . we get 



so that for 9 > 0, 



R(r)\\ c{H e^< {i+ °^ 1+e/2 , 0<9<2, 



/OO — IQ. 
e itT R(r) dr G C(R+, £(# e , L 2 )), a > 0. 
oo— ia 
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The expansion for U(t) stated in Theorem 1.2, with an error estimate (jl.8J) in C(H , L 2 ), 
is then obtained by deforming the contour of integration in the right hand side of (|2.9|) . 
exactly as in the proof of Theorem 1 of [26J . 

We now come to prove Theorem 1.2 in the Banach space C(L 2 ,L 2 ). In doing so, we 
shall make use of an argument which comes essentially from [IB], there given in the 
context of the wave equation in the exterior of a nontrapping obstacle. See also [3j. 
We shall transform the initial value problem (jl.lj) into a non-homogeneous problem. 

Let e G (0, A(oo)) be such that there are no eigenfrequencies r with the imaginary 
part Imr = A(oo) —e, and let Ti, . . . tn be the eigenfrequencies with Imr,,- < A(oo) —e. 
We introduce the spectral projector onto (&f=i E Tj , 

n =2^/< T -^ 1<iT - 

where 7 consists of the union of small positively oriented circles centered at the r^-'s. 
When x G TC, it is true that 

(r - (/ - n) x 

is a holomorphic function in a neigbourhood of the set where Imr < A(oo) — e, with 
values in Ti. We shall now derive weighted decay estimates for e ltA (I — Ii)x. When 
doing so, we notice first that Proposition 3.1 together with an interpolation argument 
shows that 

R( T ) = O(l) : L 2 -> H\ |Rer|>0(l) (3.11) 

when Imr < A(oo) — e. Let ip G C°°(R;R) be an increasing function such that 
ip(t) = 0, t < 0, and ip(t) = 1, t > 1. We consider the function ip(t)u(t,x) where 
u(t, x) solves the evolution problem (jl.lj) with the initial data (wo, u i) = (I — n)x, 
when x = (0, /), / G L 2 . Then 

(-D 2 - A + 2ia(x)D t ) i>(t)u(t, x) = g(t, x), (3.12) 

where g G C(R; L 2 ) is given by 

g(t, x) = ip'it) (2a(x)u(t, x) + iD t u(t, x)) + iD t (ijj'(t)u(t, x)) . 

The support of g(-,x) is contained in [0, 1] so that g G £ 2 (R; L 2 ), and it follows from 
the Hille-Yosida theorem that 

\\g\\mn^)<0(l)\\f\\ L 2. (3.13) 

Taking the Fourier transform in t of (j3.12j) . !Ft-* T i we S e t 

P(r)F t ^(t)u(t, x) = F^g(t, x), Ft^g{t, x) G L 2 (R; L 2 ). (3.14) 

Now an application of (j2.4j) shows that i?(r)jF^ r ^(t, x) = P(T)~ 1 J r t ^ T g(t,x) is the 
first component of 
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so that it depends holomorphically on r in the set where Imr < A(oo) — e, and takes 
values in H 1 . 

An application of Parseval's formula to (|3.14j) together with (J3.11|) and (J3.13j) gives 
that 

||e (A(oo) - £) V(t) M (t)|| L2(R; Hi) (3.15) 

— I I R( T )3 7 t-*Tg{t- l ■) |r=-+i(A(oo)-e) I \h 2 (R.;H 1 ) 

< 0(1) 1 1 ^(i, •) |r=-N(A(oo)- e ) ||l2(R;L2) < 0(1) 1 1 / | 

This gives an exponential decay estimate in the integrated form. When deriving esti- 
mates that are pointwise in time, we consider the non-homogeneous equation satisfied 
by ip T (t)u(t,x), where ip T {t) = %fj{t — T), T > 1. We have 

{-D 2 t - A + 2ia(x)D t ) i[) T (t)u(t, x) = g T (t, x), 

where gxit^x) = (V't(^) + 2a(x)ip' T (t))u(t,x) + 2iip' T (t)D t u(t,x). Standard hyperbolic 
estimates (see [^j) show that 

rT+l / fT+l \ V2 

\\u(T + l,-)\\ L 2<0(l) j \\g T (t,-)\\ L 2dt<0(l)U \\g T (t,-)\\hdtj , 



and therefore 

fT+l 



1 1 u(T + 1, •) | \ 2 L2 < 0(1) (J + (|| u(t, •) \\ 2 L2 + 1 1 D t u(t, •) | |i 2 ) . 
Now it follows from ()3.15|) that 



(3.16) 



POO 

(t,-)\\l*dt < e -2(A(oo)- e )T / e^^VWIIw^OII^^ (3-17) 

< 0(l)e" 2 ^t°°>- e ^|| / ||i 2 , 

and we shall now verify that the corresponding estimate holds when u(t, •) is replaced 
by D t u(t, ■). In doing so, we observe that ()3.14j) gives 

P{r)T t ^ T (D t (il>(jt)u(t,x))) = T^ T g{t,x), 

and since 

tR(t) = 0(1) : L 2 -> L 2 , |Rer| > 0(1), Imr G [-1, A(oo) - e], 
repeating the previous arguments we see that 

e 2(A(oo)- £)i ! ! Dmt)u ^ .))\\l 2 dt< 0(1) 1 1 f\\\ 2 . 
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Since ip' is supported in [0, 1], we get 

POO 

/ e 2 < A <°°>- e V (t) 1 1 D t u(t, -)\\hdt< 0(1)| | /Hia, 
Jo 

and therefore 

r+i 

IIA«(*,Olli^<o(i)ll/llia- ( 3 - 18 ) 

Combining (|3.16j) with (|3.17|) and ()3.18|) we conclude that 

\\u{t,-)\\ L ,<0{l)e-^-^\\f\\ L 2. (3.19) 
In order to derive an expansion of Theorem 1.2, we now only have to consider 

-i poo— ia 

e itA Ux =— e itT (r - AY 1 Uxdr, a > 0, x = (0, /), 

which is given by the sum of the residues of the integrand at the poles Tj, j — 1, . . . N. 
If tq is such a pole, then near To, we can write for some M G N, not exceeding the 
multiplicity of r , 

modulo a function holomorphic in a neighborhood of r . Here il is the spectral 
projector onto E To and Dj = (A — ro) j rio. The corresponding contribution to the sum 
of the residues is of the form e itTo times a polynomial of degree at most M in i with 
values in C(TC, D(A°°)). In particular, when r = 0, it is true that the corresponding 
generalized eigenspace Eq is spanned by (1, 0) so that the pole at is simple, and the 
spectral projector n is given by 




where 

/ 2ia(x)u (x) dx — J U\(x) dx 
J 2ia(x) dx 

Combining these remarks together with ()3.19|) . we obtain the statement of Theorem 
1.2. 

4 Proof of Theorem 1.4 

In this section we shall work under the assumption that the geodesic flow on M is 
periodic. To fix the ideas, we shall assume throughout that exp (ttH p ) (x, £) = (x,£), 
(x, £) G p~ 1 (l). Our purpose is to derive an expansion for the damped wave propagator 



3 2(A(oo)- £ )T 



7 
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U(t) on M. In doing so, we shall make use of the methods of Let us recall that 
in jT2| it was established that 

oo 

Spec (A) D {r G C; Re r > 0} C (J I k , 



k=i 



where 



4 k 4 A; 



+ i[O,0(l)], aGZ, A; = 1,2,.... 



We shall now complement this result by deriving resolvent estimates in the gaps 
between the rectangles l k . 

Proposition 4.1 There exists C > such that when Im r G [—1, 2| | a ||i°° + l] we have 



C 



\\R(r)\\c(L^) < YTM' l ReT l- C '' dist (|Rer|,U^ =1 J fc )) > 
HereJ k = [k + f - £ , A; + f + £ ] . 



1 



(4.1) 



Proof: Applying the semiclassical reduction of Section 2, we see that in order to prove 
1)4.1]) it suffices to establish the existence of C > and ho > such that for all 
h G (0, /&o] we have that 

C 



\C(L 2 ,L 2 ) 



< — . 



(4.2) 



when Imz = 0{h) and dist (Re z, U h % 1 Jk(h)) > ft,/C. Here the intervals 



Jk{h) 



# (k+^) 2 -o{h%tf (k+^y+o(h 2 ) 



k 



h 



have width 0(h 2 ) and are separated from each other by a distance which is ~ h. We 
also recall from [2Z1 and ^2] that 



Spec(P)n [a, f3] C U™ =1 J k (h). 



(4.3) 



When proving the bound 1)4.2]) . we follow the suggestion of the remark given after 
Theorem 3.2 in ^2j, which in turn is based on the methods of Sections 2 and 3 of 22j. 
Thus, repeating the arguments of the remark in [12], we find that there exists an elliptic 
selfadjoint /i-pseudodifferential operator A = Op h (e 9 ), with g = g(Rez) G S' _1 (T*M), 
such that 



A~\P + ihQ(z))A = P + ihQ + h 2 R(z), R(z) G Op h (S L 
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where Q = Op h (q) is selfadjoint, and S° 3 q is such that 

1 f T 

q=(q)T=7pJ qoexp(tH p )dt on p (Rez). 

Here q = q(Re z) = 2a\/Re z and T = T(Re z) is the common period of the closed 
ifp-trajectories in p _1 (Re z). As in the proof of Proposition 3.1, we choose g satisfying 
H P g = q — (q)r on p~~ 1 (Rez), and we may take for example, 

1 f T 

g = — J tg(exp (tH p )) dt on p (Rez). 

Since 

A, A' 1 = 0(1) : L 2 — >■ L 2 , 
it suffices to show that the bound (|4.2j) holds for the resolvent of 

V = P + ihQ + h 2 R(z). 

Now H p q = on p _1 (Re2;), and we have 

H p q = k(p — Re z), 

with k G iS 1-1 . If we let K G Op h (S'~ 1 ) have fc as a principal symbol, then 

z/j [p, Q] = h 2 K(P- Re ^) + /i 3 S, (4.4) 

where 5 G Op h (S°). Now, if A and B are bounded selfadjoint operators, we have 

\\{A + iB)u || 2 = \\Au\\ 2 +\\Bu\\ 2 + i{[A, B]u, u) 

and using this identity together with (|4.4|) we get, 

2|| (V-z)u\\ 2 > || (P + ihQ - z)u\\ 2 - 0{h 4 )\\u\\ 2 (4.5) 

> || {P-Rez)u\\ 2 + ih{[P,Q]u,u) -0{h A )\\u\\ 2 

> || (P-Rez)u|| 2 -/i|| [P,Q]u||||m||-C(/i 4 )||m|| 2 

> || {P-Rez)u\\ 2 -0{h 2 )\\ (P-Rez)wlHMI - C(/i 3 )|| m || 2 

> || (P-Re^)M|| 2 -C(/i)|| (P-Rez)M|| 2 -0(/i 3 )||M|| 2 

> ^\\(P-Rez)u\\ 2 -0(h 3 )\\u\\ 2 , 

provided that h is small enough. When z is such that Imz = 0(h) and 

dist (Rez, Uf =1 J k {h)) > h/C, 
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for some sufficiently large C > 0, it follows from ()4.3|) together with the spectral 
theorem that 

h 

1 1 (P — Re z)u || > — 1| w ||, 
and using this in ()4.5|) we obtain, for h small enough, 

2\\(V-z)u\\ 2 >^\\u\\^-0(h 3 )\\u\\^>^\\u\\ 2 . 

This immediately implies (j4.2j) with a new constant C. The proof is complete. □ 

Repeating the arguments of Section 3 and using also 1)2.70 . we find that for r as in 
Proposition 4.1, as well as for all r with Imr > 2|| a \ \oo + 1, we have 

\\R{t)\\cW)< (t+^L/a . 0<tf<2. (4.6) 

We now come to prove the propagator expansion. When G N is sufficiently large, 
we consider the positively oriented contour 7 given by the following four line segments: 
7l = { x - i ; -fc - a/4 - 1/2 < x < fc + a/4 + 1/2}, 72 = {k + a/4 + 1/2 + iy; -1 < 
2/ < 2|| a ||oo + 1}, 73 = {a; + i(2||a|| 0O + l);-A;-a/4-l/2 < x < k + a/4 +1/2}, and 
74 = {-£;-a/4-l/2 + n/;-l < y < 2| | a | ^ + 1}. Let / G H e , 9 > 0. An application 
of the residue theorem allows us to write 



! e^R(r)fdr= £ e^ r (t)/ 

TeSpcc(^), |Rcr|<fc+l/2+a/4 

/ e itT R(r)fdT- j e ltT R{r)fdT. 



'72U74 «^73 

Here p T is a polynomial in t with values in C(L 2 , H s ), for any s6R. The L 2 -norm of 
the integral along the segment 72 does not exceed 

^ , I MM O.A/9 U 4/> 



-1 



^ + |Jfe|)l+9/2 



which tends to zero as k — > 00. Similarly it follows that the contribution coming from 
the segment 74 vanishes as A; — ► 00. When treating the contribution coming from the 
segment 73, now with k = 00, we notice that in view of (J4.6)) and since 9 > 0, the 
contour 73 : Imr = 2| | a | + 1 can be replaced by Im r = C for any C > 2\\a\\ 00 + 1. 
Letting C — > 00, we conclude that the corresponding contribution to the integral 
vanishes. We get with k ^> 1, 



00 



U(t)f= e Ur p T (t) + J2 U k(t)f, whent>0, / G H e } (4.7) 

TeSpec(yl), |ReT|<fc +o/4 k=k 
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where 

U k (t)= J2 e "X^)+ J2 e itT Pr,k(t), 

T£Spec(A),T<=I k reSpcc(^), -rel k 

and p Tt k{t) is a polynomial in t with values in C(L 2 , H s ), for any sGR. An application 
of ()4.6|) together with the residue theorem shows that for each fixed t we have 

||i/ fc (t)|| jC(HfliLa) = 0(l)e* 1 



Jfel 



1+0/2 ' 



The absolute convergence of the series (|4.7jl in C(H 9 ,L 2 ), 9 > 0, follows, and this 
concludes the proof of Theorem 1.4. 



5 Proof of Theorem 1.3 

As in Section 1, we consider the operator 

V = P + ihQ(z), P = -h 2 A, Q(z) = 2a(x)y/Z, z G [a, /3] + i[- 7 , 7], 

for 0<a<l</3, 7>0. In the following discussion, we shall assume that the if p -flow 
possesses a simple closed trajectory 7 : [0, T] — ► p _1 (l), which is nondegenerate elliptic 
and such that 

7r( 7 ([O,T]))nsupp(a) = 0. (5.1) 

Here tc : T*M — > M is the natural projection. Moreover, we shall assume that the 
linearized Poincare map P 1 of 7 is A^-elementary, for some > 4. We shall establish 
the existence of a sequence of eigenvalues of the operator A converging rapidly to the 
real axis. 

Under these assumptions on 7, it is proved in [H] that there exist quasimodes for 
P = P(h) associated with 7 and having polynomially small errors. To be precise, it 
follows from the results of [0] that there exists a set H C (0, 1] with e H, the closure 
of H, and < ao < a(/i) < b(h) < &o < 00, two functions on H, such that for any 
h G H there exists m(h) G N, Ej(h) G [a(h),b(h)], and Uj(h) G L 2 , j = 1, . . .m(h), 
such that 

(P(h) - EjQi)) Uj (h) = R (h) in L 2 , j — 1,... , m{h), 

and (uj(h)\uk(h)) = 5jk + Ro(h), where Ro(h) = 0(h°°). Furthermore, it follows from 
the construction of that if A G Op h (S°) is such that WF h (A) is contained in a small 
neighborhood of 7 ([0, T]) and A = I microlocally near 7([0,T]) then 

(/ - A)uj(h) = 0{h°°) in L 2 , j = 1, . . . m{h). 

From the assumption (|5.1jl we infer therefore that 

(V - Ejih^Ujih) = R(h) in L 2 , j = 1, . . . m(7t), 
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where R(h) = 0(h°°). 

Remark. The construction of Q\ gives quasimodes of P associated with a family of KAM 
tori of the Poincare mapping, near the closed trajectory 7. For analytic manifolds, the 
construction of |E] produces quasimodes with exponentially small errors. We may also 
recall that under the stronger assumption of the iV-elementarity of P 1 for every JVgN, 
the existence of quasimodes with polynomially small errors can be established also by 
the method of Gaussian beams — see and the references given there. While the 
construction of |E] gives quasimodes of positive mass, given by the volume of the union 
of the KAM tori, it is known that the quasimodes obtained by the method of Gaussian 
beams have mass zero — see [7j. 

For future reference, we notice that for z G [a, (3] + i[— 7, 7], we have Im (Vu\u) > 0, 
and therefore 

\\{V- z)- 1 |U (x2)i2) < — lmz<0. (5.2) 

\Lmz\ 

The following result implies Theorem 1.3. 

Theorem 5.1 For any positive function S(h) with De~ D l h < S(h) = 0(h°°), for some 
D > 0, and such that S(h) ^> R{h), and every k e N, there exists h(S,k) > such 
that for H 3 h < h(S, k), the spectrum ofV = V(h) intersects the set 

[E j {h)-Qh k ,E J {h) + Qh k }+i[0,2S{h)h- n - 1 ], j = 1, . . . m{h). (5.3) 

When proving Theorem 5.1, we shall start by establishing the following result, 
which gives an a priori exponential estimate of the resolvent of V, in a set of size ~ h, 
outside the union of small neigbourhoods of the eigenvalues. 

Proposition 5.2 Let E 6 [a + 1/C , j3 — 1/C ] for some C > 1. For every L > 
there exists A > and h > such that for < h < h we have 



||(P-,)-'|U ( ^,<^exp(^-»log^) 
for \z-E\< Lh, dist (z, Spec('P)) > g(h). Here < g(ti) < h. 



(5.4) 



Proof: We shall use the idea of [I7j and [221 °f exploiting a finite rank perturbation to 
create a gap in the spectrum of P around E. Recall from |22] that for any C > there 
exists a selfadjoint operator P = P + K, with the same domain as P, such that 

[E -Ch,E + Ch] n Spec (?) = 0, 

and 

\\K\\ c(L2 , L2) <Ch, \\K\\ tr <C(C)h 2 - n . 
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Here || • || tr denotes the trace class norm. We then introduce 

T = P + ihQ{z), 

and writing 



we see that 



V - z = (P - z)(I + ih(P - z)- l Q), 
{V-z)-^^- in C{L\L% 



(5.5) 



when z is such that |Re;? — E\ < Lh, provided that C is large enough. Now we write 

V-z= (V-z)(I- (V-z)- x K), (5.6) 

and let us introduce 

D{z) = det(J- (V-z)- l K), 

which is a holomorphic function in the open disc D(E, Lh) = {z G C; \z — E\ < Lh}, 
whose zeros are precisely the eigenvalues of V in this set. Since the trace class norm 
of (V - z)~ 1 K is 0(l)h l - n } we get 

\D(z)\ < exp(C(l)/i 1 - n ), z G D(E, Lh). (5.7) 

An application of Theorem 5.1 from Chapter 5 of JT] shows that when |Re z — E\ < Lh, 



(I-(V-z)- L K)- l \\ c{L 2 jL2) < 



det(/+ (V-z)-^ 



D z) 



< 



\D(z)\ ' 



(5.* 



and when deriving bounds on the resolvent of V, we shall therefore have to estimate 
|.D(z)| from below, away from its zeros. In doing so we remark that when \z — E\ < Lh 
and |Im^| > 0(l)h then 



(V-z) 

and from (|5.6|) we get that in this set, 



h 



-l 



I-(V- z)- L K =(V- z)-\V -z) = I+(V- z)- L K 



with 

It follows that 



(v- z y l K\\ tT <o{i)h 



l-n 



\D(z)\>e-°^ hl - n . 



z G D(E, Lh), \lmz\ > 0(l)h. 



(5.9) 
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Let now Z\,...z^ be the zeros of D(z) in D(E, Lh). From j22j we recall that iV 
£>(l)/i 1 - n . We factorize 

D(z) = G(z)B{z), z G D(E, Lh), 
where G is holomorphic and non- vanishing in D(E,Lh), and 

N 

B(z) = J[b Ss (z) 



is the Blaschke product on D(E, Lh), so that 

b z Az) =b Zj 



z — E\ z — Zj 



Lh ) Lh(l-(^)(^ E 



Lh 

Here 

K{z) = - — — , \z\ < 1, \w\ < 1, 
1 — wz 

stands for the standard Blaschke factor for the unit disc. An application of basic 
factorization theorems in the Banach algebra of bounded holomorphic functions in a 
disc, or simply using the maximum principle together with ([5.7)1 . allows us to conclude 
that 

\G(z)\ <e° (1)hl ~ n in D(E,Lh). (5.10) 
Taking (|5.9|) into account we find that for |Im^| > 0(l)h, 

\G(z)\>e-°^ hl - n , 

An application of the Harnack inequality to the harmonic function log|G(z)|, in the 
form given in Chapter 1 of ^7j, shows that 

|log|G(z)|| < 0(l)h}- n , in D(E, Lh), 

after decreasing L slightly. It is easy to see that 

\B{z)\ > e -° {1)hl ~ niog ^), z G D(E, Lh), dist(z,Spec(V))>g(h), 

and therefore for such z's we get, 

|^(^)|> e -° (1) ^" los ^. (5.11) 

Combining ()5.6|) and ()5.5|) with (|5.8|) and ()5.11|) . we complete the proof of ([5.4)1 . □ 

Remark. The resolvent bound of Proposition 5.2 together with its proof is closely 
related to the abstract resolvent estimates for weak non-self adjoint perturbations of 
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selfadjoint operators, developed in Chapter 1 of [T7j. See also j^j and jTUj for further 
references and applications of similar bounds in the theory of resonances and semiclas- 
sical pseudospectral theory. 

Remark. For the following argument, it would have been sufficient to derive an expo- 
nential bound on (V — z)^ 1 in a bounded h- independent set, rather than in a set of size 
~ h. Indeed, if Q CC {z G C; Rez > 0} is a bounded open /i-independent neighbour- 
hood of some E > 0, then a straightforward adaptation of the proof of Proposition 5.2 
shows that for all h small enough and z G Q, we have 

II (V-z)- 1 \\cw) < 0(l)exp(0(l)h- n log(l/g(h))), dist (z, Spec(P)) > g(h) > 0. 

(5.12) 

See also Lemma 6.1 of 10 . The exponent h n in ()5.12j) reflects the fact that the 
number of eigenvalues of V in Q is 0(h~ n ) — see [22] for the much more precise result 
giving the Weyl asymptotics. It seemed to us nevertheless worthwhile, and indeed, 
of independent interest, to obtain a sharper result in the smaller set, whose size is 
dictated by the strength of the non-selfadjoint perturbation. Notice also that in the 
one-dimensional case, n — 1, Proposition 5.2 gives a polynomial bound on the resolvent 
of V. It follows that in this case, only eigenvalues can produce a super-polynomial 
growth of the resolvent. 

We now come to prove Theorem 5.1. When doing so, we follow the argument of j2H| 
and argue by contradiction. Theorem 5.1 is then obtained by applying the argument 
of as it stands, since an application of the semiclassical maximum principle, Lemma 
2 of [2H], to (V — z)~ x is now legitimate, in view of Proposition 5.2, (|5.12|) . and the 
bound (|5.2|) . The proof of Theorem 5.1 is complete. 

Remark. Following j2B], it is possible to sharpen the statement of Theorem 5.1 and to 
estimate the number of eigenvalues of V in the set 

[a(h) -0{l)h k ,b{h) +0{l)h k ] +i[0,2S(h)h- n ~ 1 }, keN. (5.13) 

Below we shall indicate how to do so, by adapting the argument of [2S1 to the present 
setting. To this end, let us introduce the operator 

= ( —h 2 A 2 Lk ) = (o h) hA {l I)'* ■■H 1 xL'^H 1 xL'. (5.14) 

It is then true that z G Spec(P) if and only if A = \fz G Spec(A(h)). Using the 
factorization in (|5.14|) together with Proposition 2.1, we verify that 

IK^-A)- 1 !!^)^^^, ImA<0. 

Here 7i = H 1 x L 2 , and we equip H 1 with the corresponding /i-dependent norm. It is 
also easy to see, using Proposition 5.2, that the resolvent of A(h) enjoys an exponential 
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type bound of the same form as (V — z) l , when away from the eigenvalues. If we now 
use the quasimodes for P and set 

U ' [h) = ( Ef^uAh) ) £ *' > = L -•»(*). < 5 ' 15 » 

then a simple computation shows that 

(A(h) - E) ,2 {h))Uj{h) = R(h) in H, 

and {Uj(h)\U k (ti)) H = (1 + Ej(h) + E k (h))5 jk + R(h), where R(h) = 0{h°°). Here we 
have used that the modes Uj(h) are orthonormal in L 2 and ortogonal in H 1 , modulo 
an error which is 0(h°°). 

An inspection of [23J combined with the preceding remarks shows that the counting 
argument of Theorem 1 of j2H| applies to the operator-valued function h(A(h) — A) -1 
without any change, and demonstrates that for any positive function S(h) such that 
max(De~ D / h , h~ n ~ 2 R(h)) < S(h) = 0(h°°), D > 0, and any k e N, the rectangle 

[a{h) 1/2 - 6h k , b{h) + 6h k ] + i[0, 2S(h)h- n - 1 ] 

contains > m(h) eigenvalues of A(h), for h small enough. It follows that the set (J5.13|) 
contains at least m{h) eigenvalues of V . 
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